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A POPULATION MODEL FOR FLEAS OF THE 
GRAY-TAILED VOLE, MICROTUS CANICAUDUS 

MILLER' 



Richard G. Robbins^, G. David Faulkenberry^ 

ABSTRACT: Negative binomial probability distributions are successfully fit to observed 
frequency distributions of fleas on the gray-tailed vole, Microtus canicaudus Miller, in 
western Oregon. The host flea population is thus shown to be contagiously dispersed. 



This is the second in a projected series of papers on the population 
dynamics of the fleas that parasitize the gray-tailed vole, Microtus 
canicaudus Miller, a field mouse that occurs in grassy, uncultivated areas 
between the Cascade and Coast Ranges of western Oregon and Washington. 
Earlier (Faulkenberry and Robbins, 1980), the odds ratio was used to 
measure the degree of association between different flea species on this 
vole. However, ecologists also recognize the predictive value intrinsic to a 
mathematical description or model of the dispersion of organisms in space, 
and in medical entomology such spatial relationships between hosts and 
their parasites are often expressed by the negative binomial (or Pascal) 
distribution (Bliss, 1953; Crofton, 1971a, b). The negative binomial is a 
theoretical frequency distribution or probability distribution that may be fit 
to frequency distributions observed in nature . Thus, Fisher (1941) success- 
fully fit the negative binomial to Milne’s ( 1 943) frequency distributions for 
the tick Ixodes ricinus (Linnaeus) on sheep in England, thereby mathe- 
matically describing the dispersion of the parasite population on a 
particular host. More recently, Easton (1973) and Easton and Goulding 
( 1 974) fit the negative binomial to data for Ixodes angustus Neumann and 
larval trombiculid mites (chiggers) on various hosts in western Oregon. Yet, 
with the possible exception of Williams’ (1964) work on J.L. Harrison’s 
data for Xenopsylla astia Rothschild and X. cheopis (Rothschild) on the 
murine rat Bandicota bengalensis (Gray) in Rangoon, Burma, the present 
paper is the first published record wherein the negative binomial is 
successfully fit to data for fleas. 
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Application of the Negative Binomial Distribution 



Specifically, the negative binomial can arise ( 1 ) as a result of a series of 
exposures to parasites in which each exposure is random and the chances of 
acquiring parasites differ at each exposure, (2) as a result of nonrandom 
dispersion of parasites, (3) as a result of differences between individual 
hosts that affect the chances of acquiring parasites, and (4) as a result of a 
change with time in the chances of being parasitized (Crofton, 1971a). The 
negative binomial is described by two parameters, the mean and the 
exponent k, and is fit to contagious (over-dispersed or clumped) populations 
in which the variance always exceeds the mean. With combinatorial 
notation, the distribution takes the form: 



p(x) = 



k + x - 1 
k - 1 



(l-0) x (9 k 



x =0, 1, 2, . . . 



where the mean is k0/(l-0) and the variance is k0/( 1-0)2. When the 
dispersion parameter k is small, the degree of clumping is great and the 
variance is much larger than the mean, but as k increases the distribution 
approaches a Poisson series in which the variance equals the mean. 

Among the many methods of estimating k from sample data are the 
following three given by Southwood (1978): 
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mean number of parasites per sample (host) 
variance 

total number of samples 

number of samples bearing no parasites 

Napierian logs 

the sum of all frequencies of sampling units bearing more than x 
parasites 



Unless the mean is quite low, formula ( 1 ) is not reliable when k < 3, that is, 
when populations show a moderate degree of clumping. For this reason, 
formula (1) is usually used only to describe low density populations. 
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Formula ( 2) is applicable to most populations with very small means but to 
large ones only when there is extensive clumping. Generally, about one- 
third of the hosts must be free of parasites if the mean is below ten, and as the 
mean increases greater numbers of hosts must fall into the zero class. 
Formula (3) is considered a better estimator than either ( 1 ) or (2) but is also 
subject to bias when the mean is very small and k very large. 

Once values of k have been obtained, the negative binomial population 
model can be fit in three ways to the observed frequency distribution. First 6 
is estimated by k/(k + x), then expected frequencies are calculated by: 



where P x is the probability of a host bearing x parasites and the values x! 
and y(k) are obtained from tables of factorials and of log gamma functions 
respectively (Southwood, 1978). Expected and observed frequency distri- 
butions are compared by a chi-square which has three fewer degrees of 
freedom than the number of comparisons that are made. Chi-square 
comparisons are considered strongest when expected frequencies are at 
least 1. 



Table 1 compares the k, chi-square, and degree of freedom values for 
negative binomial distributions fit to observed frequency distributions of 
Atyphloceras mnltidentatus (C. Fox), Catallagia charlottensis (Baker), 
and fleas collectively on 377 comparable voles taken from ecologically 
similar sites to the north, south, and west of Corvallis, Oregon (for an 
account of sampling techniques and a list of all flea species obtained during 
this study, see Faulkenberry and Robbins, 1980, or Robbins, 1976). 
Because formulas (2) and (3) for the calculation of k must be made to 
balance by the method of iteration, all models were generated by computer. 
In order to obtain the best possible fit, each of the three observed frequency 
distributions was assigned two values of G, the expected frequency below 
which all classes were pooled. Of course, in the case at hand, pooling was 
only possible when G was set at 1 because all classes had expected 
frequencies greater than 0. In addition, each distribution was subjected 
once to truncation by removal of the zero class, which is by far the largest 
class in every case. In Table 1 the best chi-square value in each category is 
underscored. 

It is clear that the fleas of the gray-tailed vole are very contagiously 
dispersed. All three formulas for k led to the best fit an equal number of 
times and only once was the difference between expected and observed 



y (k + x) 
x! y(k) 
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frequencies significant at P < 0.05 . On the other hand, the fit of the negative 
binomial to the frequencies for zero and one flea per host was seldom close 
— an indication of top-heaviness in these classes. Consequently, the best 
distributions were almost always truncated. 



Table 1 . Comparison of k, chi-square, and degree of freedom values for all negative binomial 
distributions. 



G = 1; 
not truncated 



G =0; 
truncated 



G =0; 
not truncated 



Atyphloceras 


kl 


= 


0.41384 


kl 


= 


1.50700 


kl 


= 


0.41384 


multidentatus 


X2 


= 


2.1868355 


X2 


= 


1.5863500 


X2 


= 


2.8879404 




df 


= 


3 


df 


= 


7 


df 


= 


8 




k2 


= 


0.35500 


k2 


= 


2.08000 


k2 


= 


0.35500 




X2 


= 


1.8917037 


X2 


= 


1.9245513 


X2 


= 


2.5099838 




df 


= 


3 


df 


= 


7 


df 


= 


8 




k3 


= 


0.36202 


k3 


= 


1.74250 


k3 


= 


0.36202 




X2 


= 


1.9010235 


X2 


= 


1.6964298 


X2 


= 


2.5288784 




df 


= 


3 


df 


= 


7 


df 


= 


8 


Catallagia 


kl 


= 


0.36196 


kl 


= 


0.46719 


kl 


= 


0.36196 


charlottensis 


X2 


= 


12.3286353 


X2 


= 


6.4885838 


X2 


= 


12.5538454 




df 


= 


6 


df 


= 


7 


df 


= 


8 




k2 


= 


0.45250 


k2 


= 


0.40000 


k2 


= 


0.45250 




X2 


= 


12.3763963 


X2 


= 


6.0815673 


X2 


= 


12.7556870 




df 


= 


6 


df 


= 


7 


df 


= 


8 




k3 


— 


0.42405 


k3 


= 


0.40788 


k3 


= 


0.42405 




X2 


= 


11.8824064 


X2 


= 


6.0885154 


X2 


= 


12.1912021 




df 


= 


6 


df 


= 


7 


df 


= 


8 


All species 


kl 


= 


0.50407 


kl 


= 


0.52748 


kl 


= 


0.50407 




X2 


= 


14.7781608 


X2 


= 


15.6730234 


X2 


= 


20.7273663 




df 


= 


10 


df 


= 


17 


df 


= 


18 




k2 


= 


0.67375 


k2 


= 


0.70375 


k2 


= 


0.67375 




X2 


= 


12.0033085 


X2 


= 


21.4046443 


X2 


= 


32.1796904* 




df 


= 


9 


df 


= 


17 


df 


= 


18 




k3 


= 


0.62230 


k3 


= 


0.62526 


k3 


= 


0.62230 




X2 


= 


11.6262368 


X2 


= 


17.2498995 


X2 


= 


25.7152032 




df 


= 


9 


df 


= 


17 


df 


= 


18 



*Significant at P < 0.05. 
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